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IRAS LRS spectra (IRAS team, 1986) of M Mira variables have shown a large variation in the 
appearance of the 9.7 |im silicate feature, which is correlated with the shape of light curve (c/., Vardya et 
al., 1986). Onaka and de Jong (1987) and Onaka et al. (1988) have studied the LRS spectra of about 100 
Mira variables by using simple dust shell models containing mixtures of silicate and aluminum oxide dust 
grains. They have shown that the aluminum oxide grains account for the observed broad feature around 12 
pm and that the variation of the spectra can be interpreted in terms of the variation of the temperature at the 
inner boundary of silicate dust shell. They have proposed silicate mantle growth on aluminum oxide grains 
as a possible explanation for the results. In this report, we calculate model spectra taking account of sili- 
cate mantle growth and investigate the physical parameters which may determine the appearance of the 9.7 
pm feature in M Mira variables. 

In the model calculation it is assumed that aluminum oxide grains are already formed at the bottom of 
the circumstellar envelope because of their high condensation temperature (- 1500 K). The growth of sili- 
cate mantle and the motion of gas and grains from r =r 0 , where the mantle growth starts, are investigated. 
Sticking and sputtering processes due to the relative motion of grain to the ambient gas are taken into 
account. The thermal velocity is assumed to be negligible to the drift velocity. Acceleration by radiation 
pressure is considered in the gas motion equation. The formal solution is integrated to obtain the emergent 
spectra. Physical conditions inside r 0 are regarded as boundary conditions. Observed spectra are com- 
pared to model spectra to investigate the conditions at the bottom of circumstellar envelope. In modelling 
the envelope, a parameter C/ is introduced to take account of the density fluctuation of the envelope 
phenomenologically. 

Silicate mande growth on an aluminum oxide grain can be written as 
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with A=ala 0 , R=r<Jr, U=v/v 0 , V d =v/v d , and A f =a f la 0 . Here, a is the grain radius including core and 
mande, a 0 the core radius, r the distance from the star, v the gas velocity, Vq=v at r=r Q , and v d the drift 
velocity of the grains relative to the ambient gas. The parameter aj is determined by the abundance rado 
of silicate to aluminum oxide, being about 3. The last term in the parentheses represents the erosion due to 
sputtering: f - t is the relauve cosmic abundance of element i and S t the sputtering yield taken from Kwok 
(1975). Sputtering is negligible unless v d exceeds 20 km/sec. The parameter F„ describes the rate of man- 
tle growth and is given by 
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where M is the mass loss rate, / the fracuon of the condensible material (- 310 ), Q. the volume of the 
monomer (~ 5 TO -23 cm 3 ), p the molecular weight (- 1.4), and m H the hydrogen mass. 

The change of the velocity due to the radiation pressure may be given by 
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where Q p is the radiation pressure coefficient of the grains and the parameters g and F v are defined as 
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Here L is the stellar luminosity (~ 2.5- 10 37 erg/sec), M the stellar mass (~ 210 33 g ),f g the mass fraction of 
aluminum oxide (~ 9.6- 10 5 ), d the density of aluminum oxide (~ 4 g/cm 3 ), c the light velocity, and G the 
gravity constant. 

The drift velocity variable V d is given, if the thermal velocity is much smaller than v d , by 
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The emergent flux is obtained by the integration of the formal solution (Chandrasekhar, 1960): 
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Here T. is the stellar temperature (~ 2500 K), R. the stellar radius (~ 310 13 cm), D the distance to the star, 
and x x the extinction optical depth. The input parameter p is normalized by R, . The dust temperature 
T d {r) is calculated from the energy balance: 
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and the volume emission coefficient A (r ) is given by 
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where r x is the inner boundary of the aluminum oxide dust shell (~ 3.8- 10 13 cm) and the suffix o represents 
the value at r=r 0 . The optical constants for aluminum oxide are taken from Eriksson et al. (1981) and 
those for silicate referred to Day (1979). 

In Figure 1, the radial variation of A , V, v d , and T d are shown for models with different a 0 and M . 
As seen in Figure 1, the changes of V and v d are small. The gas velocity decreases slightly at the beginning 
in order to adjust the boundary conditions (v=v 0 ) and the decrease is probably not realistic. However, it is 
small and does not make significant effects to calculated spectra. The temperature profile weakly depends 
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on the model parameters. Thus, calculated spectra depends mostly on the parameter F R and lc 0 . If the 
grain size is much smaller than the wavelength in question, Q p and Q abs are approximately proportinal to 
the size. Therefore, F R is proportional to M la 0 and k(r) to M: a 0 and M are the major parameters in this 
model. 
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Figure 1. a) Radial distribution of grain radius normalized by core radius, b) Normalized velocity, c) Dust 
temperature (in K). d) Drift velocity normalized by initial gas velocity. All models are calculated 
with v 0 =10 km/sec, C t = 1, and r o =6T0 14 cm. Model numbers are indicated on each curve; 1: 00=10 
nm, M~\0~ 5 M^/yr, 2: 0 o =lOO nm, M~ 10” 5 M^/yr, 3: 0 q= 2O nm, A/=10"^ Mo/yr, and 4: 0 o =lOO 
nm, M- 10~ 6 Me/yr. The abscissa is the normalized distance (r/r 0 ). 


A model grid was constructed for various sets of a 0 and Af. The best fit model parameters were 
obtained for each observed spectrum. The samples of LRS spectra are the same as in Onaka et al. (1988). 
The parameter r 0 was set to be 6T0 14 cm since it was found to give the best fit of models to most observed 
spectra. Effects of v 0 and C t were also examined. Examples of fitted spectra with LRS spectra are shown 
in Figure 2. Most LRS spectra of M Mira variables in the present sample are reproduced satisfactorily by 
the models with parameter range KT 6 < M < 10” 5 M^yr and 5 < a 0 < 400 nm). However, in some cases, 
e.g. t RR Aql, a much larger F R is reqiured to reproduce a strong silicate feature. Unless we assume a very 
fine core radius (~ 0.5 nm) it is difficult to have a good fit. 
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Figure 2. Examples of model spectra together with LRS spectra of M-type Mira variables. They are 
shown in the order of light curve asymmetry index / (cf. Vardya et ai, 1986). Model parameters 
are: (a„ (nm), M (KT 6 M©/yr), v 0 (km/sec), and C ,) = (5, 2, 10, 10) for RR Aql, (5, 2, 7, 1) for X 
Cen, (10, 2, 7, 1) for RR Sgr, (200, 10, 7, 1) for RR Per, (100, 10, 10, 1) for RS Lib, and (400, 10, 7, 
5) for R Car. 
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Figure 3. Distributions of model parameters Log( M) (M in M^yr) and Log( a 0 / C x ) ( a 0 in nm) for stars 
with different light curve asymmetry index / . 


Thus, we introduce a dumpiness parameter C x as described in equations (2) and (6). We assume that 
in these stars a density fluctuation occurs and accelerates the mantle growth. It also affects the drag force 
(equation (6)). However, owing to the optically-thin nature of the present models, effects to the transfer 
equation (7) can be neglected. If we introduce C,, the model lit becomes much improved ( e.g ., RR Aql in 
Figure 2). As inferred from equations (2) and (6) it is difficult to separate the effects of C x from those of a 0 
a priori. The change of v 0 also improves the model fit. However, its effects are small since observed ter- 
minal velocities are in a small range. 

The distribution of model parameters a 0 divided by C x and M are shown in Figure 3 for five groups 
of stars with different light curve asymmetry index / . It is clear that for stars with small / index (asym- 
metric light curve) a small core radius or a large clump is required. For models of stars with large/ index 
(symmetric light curve), somewhat larger M is obtained compared to stars with small / index. This is con- 
sistent with the observed trend of the 9.7 |im feature with / : Only stars with / <0.43 show the 9.7 |im 
feature and those with / >0.43 do not. If a 0 is small or C t is large, F R becomes large and silicate mantle 
grows quickly; emergent spectra show the silicate feature clearly. If a 0 is large, on the other hand, mantle 
growth is suppressed and the aluminum oxide feature at 12 Jim is observed. 
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According to the present model, the observed variation in appearance of the 9.7|im feature is ascribed 
to the variation of grain size of core aluminum oxide grains or to the degree of density fluctuation in the 
circumstellar envelope. Some degree of density fluctuation is necessary unless very fine grains are 
assumed, although in the present analysis it is difficult to separate the effects of C, from those of a 0 . The 
size of aluminum oxide grains is considered to be determined by the nucleation process at the bottom of the 
circumstellar envelope. According to Yamamoto and Hasegawa (1977) and Draine and Salpeter (1977) the 
particle size in the homogeneous nucleation process is determined by the ratio of the cooling time scale of 
the system to the collision time scale of the condensible gas particles. In stars with small/ index the shock 
propagates strongly and the cooling occurs rapidly, thus, small aluminum oxide grains may be formed. It 
may be also possible that further mutual collisions occur frequently in circumstellar envelopes of these 
stars and that grains are broken into small pieces (Biermann and Harwit 1980). Strong shocks may also 
produce a large density fluctuation. The observed variation of LRS spectra, thus, can be ascribed to the 
difference of the physical conditions at the bottom of the circumstellar envelope and should provide impor- 
tant information on the acceleration mechanism of mass loss process. 

The increment at 18(j.m for some stars (e.g., R Car in figure 2) may not be well reproduced by the 
present model. The dependence of silicate band strength on the temperature, which has been indicated 
experimentally by Day (1976) and suggested from the dust shell model analysis by Bedijn (1987), may 
have to be taken into account. 
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